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METHOD AND SYSTEM FOR PREDICTIVE CONTROL OF LIVE 

STREAMING VIDEO MEDIA 

BACKGROUND OF THE INVENTION 
The present invention relates generally to digital data processing of video 

information. More particularly, the invention provides techniques, including methods and 

systems, for encoding video information. 

The video information is usually transmitted from a video provider to a video 

user over a transmission system, such as a transmission line or a wireless network. The video 

information includes motion pictures, live news broadcasts, and real-time information 

i 

| updates. The video provider can be, among others, an Internet service provider, a 
j videoconference participant, an entertainment service provider, or a news agency. The video 
user is, among others, an Internet service subscriber, a videoconference participant, a video- 
enabled cellular phone user, a CODE C user, a PALM user, or a user of any other types of 
handheld devices that possess the video capability. 

The transmitted video information usually comprises thousands of frames, 
each of which is composed of numerous microblocks. Each microblock usually represents 
one pixel of the image illustrated in a frame, and comprises multiple bits of data. During the 
transmission, each frame lasts from about 25 to 30 msec in order to generate a continuous 
* motion effect for the video viewer. The unencoded video information hence needs an 

average transmission rate of millions of bits per second. In addition, because various frames 
may contain drastically different amount of data, the required transmission rate could vary 
between different frames. 

In contrast, at the hardware level, the transmission rate of the transmission 
system is usually much lower than millions of bits per second, and it usually does not vary 
with time. This limited transmission rate cannot typically support the unencoded 
transmission of the video information; therefore the video data must be encoded in order to 
reduce the required transmission rate. 

The required transmission rate of the encoded data on the average needs to be 
equal to or smaller than the transmission rate of the transmission system. However, if the 
video date are overly compressed, video images that the viewer sees after the transmission 
may not possess a satisfactory quality. Furthermore, if the required transmission rate for the 
encoded video data drops below the transmission rate of the transmission system, the unused 



transmission capacity is filled with meaningless information, and is wasted. Hence, the 
optimization of encoding parameters is critical for transmitting high quality video images 
with the limited transmission rate. 

The search for the optimal encoding parameters is complicated by the 
significant variation in data size of various frames. With current encoding schemes, encoded 
video data have file sizes that are not strictly proportional to the number of microblocks in 
each frame, but rather proportional to the spatial complexity and motion activity of the image. 
The various frames, if encoded with the same set of parameters, may require different 
transmission rates, but the transmission rate of the transmission system is fixed. Therefore 
the optimal encoding parameters vary frame by frame. 

One method for determining a set of encoding parameters for each frame 
includes a buffer that stores the encoded data prior to their transmission over the transmission 
I system. When the incoming rate of the encoded data is higher than the transmission rate of 
the transmission system, the buffer stores the incoming data. As the buffer fills, the encoding 
parameters are adjusted so that the transmission rate for the subsequent encoded data is 

j reduced. If the encoding parameters are not adjusted, the buffer will eventually overflow, 
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| and consequently the incoming data would be irretrievably lost. On the other hand, when the 
incoming rate of the encoded data falls below the transmission rate of the transmission 
system, the buffer maintains its output rate at the transmission rate of the transmission 
j system. Similarly, the encoding parameters are adjusted to raise the transmission rate of the 
j future subsequent encoded data. If the encoding parameters are not adjusted, the filled 
percentage of the buffer would drop below a predetermined level, and would have to be filled 
with bits of random data to prevent the crash of the encoding system. In short, the adjustment 
triggered by the buffer can compensate the adverse effects resulting from buffer underrun and 
overrun. 

A drawback to this method is that the buffer cannot completely solve the 
problem associated with the lack of optimization. The adjustment of the encoding parameters 
based on the buffer fullness is predicated upon the dubious assumption that future frames 
have similar complexity as previous frames. In reality, this inference of future behavior 
based on recent historical data is often incorrect. Under this assumption, a simple image 
following a complicated frame would be overly compressed at the expense of image quality, 
and a complex image following a simple frame can be under compressed and thus fills up the 
buffer. 



Another drawback is that large buffers are required. In particular, a very large 
buffer may be used to allow a significant rate mismatch between the required encoded data 
transmission rate and the transmission rate of the transmission system over an extended 
period without using a complicated rate control algorithm. However, the large buffer 
increases the cost of the transmission system, and introduces a great deal of testing. Such 
added latencies delay are unsuitable, especially for communications between parties 
responding to each other interactively, such as for a video conferencing system. 

Some other methods for solving the rate mismatch problem include making 
adjustments based on historical data and the fullness of the buffer; however, this solution 
cannot solve the problems related to the lack of intra-frame optimization. For example, under 
these methods, the modifications to the encoding parameters often create a cyclic effect 
undesirable for videoconferencing. In a videoconference, the scene is usually composed of a 
very simple background at the top of the frame, talking heads with increased motion in the 
center, and hands manipulating objects on a table top at the bottom of the frame. As an 
encoder processes the simple background at the top of the frame, the compression becomes 
more and more compact in order to compensate for the accumulated buffer fullness. By the 
time the table top is reached, the objects to which the speaker is drawing the viewer's 
attention are overly compressed, and are poorly represented. 

In order to solve the problems related to the lack of the optimization, past 
researchers have suggested various methods, such as MPEG-4 Standard (See IISO/IEC FDIS 
14496-2, Annex L), for analyzing the date content of the video frames before they are 
encoded. However, these methods require complex mathematical operations, and are not 
suitable for fast encoding on a low-speed microprocessor. 

In summary, highly desirable is a predictive encoding technique that, based on 
the complexity of the future frame, adjusts the encoding parameters and utilizes only simple 
mathematical manipulations. 

SUMMARY OF THE INVENTION 
According to the present invention, a technique for digital video processing is 
provided. In exemplary embodiments, the present invention provides methods and systems 
for processing digital video data in a predictive manner for live video streaming applications. 

In one embodiment, a method for processing digital video signals for live video 
applications may include providing video data that comprise a plurality of frames, identifying 
a first frame and a second frame in the frame sequences, and processing the information of 



the first frame and the information of the second frame to determine a quantization step for 
the second frame. The step of processing the information of the first frame and the 
information of the second frame may include calculating an average of a sum of absolute 
differences (sigmaSAD) value for the second frame, calculating a divisor value for the second 
frame, and calculating the quantization step for the second frame. 

In another embodiment, a method for processing digital video signals for live video 
application may include deciding whether to encode a frame in the I-mode before any P- 
frame encoding is accomplished. 

In yet another embodiment, a method for processing digital video signals for live 
video application may include adjusting the complexity integer so that the size of encoded 
data for anyone of the plurality of frames has approximately an equal size. 

In yet another embodiment, a method for processing digital video signals for live 
video application may include determining the locations of I-frames in the step of providing 
video data, extending frames immediately preceding the I-frames for one additional frame 
time, and skipping frames immediately following the I-frames. 

In yet another embodiment, a system for processing digital video signals for live 
video applications comprises a memory unit within which a computer program is stored. 
Functions of the computer program may include instructing the system to provide video data 
comprising a plurality of frames, instructing the system to identify a first frame and a second 
frame in the frame sequences, and instructing the system to process the information of the 
first frame and the information of the second frame to determine a quantization step for the 
second frame. 

In yet another embodiment, a system for processing digital video signals for live 
video applications is disclosed. The system may include a video providing subsystem that 
provides video data comprising a plurality of frames, a sigmaSAD calculation subsystem that 
calculates a value of sigmaSAD, and a divisor calculation subsystem that calculates a value 
of divisor. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 illustrates a block diagram for a method for processing digital video 

signals according to an embodiment of the present invention; 

Fig. 2 illustrates a block diagram for a method for processing digital video 

with the capability of advance determination of I-mode frames according to an embodiment 

of the present invention; 



Fig. 3 illustrates a block diagram for a method for processing digital video 
with the capability of ensuring an approximately same size for encoded data of various 
frames according to an embodiment of the present invention; 

Fig. 4 illustrates a block diagram for a method for processing digital video 
with the capability of maintaining high video quality while skipping frames according to an 
embodiment of the present invention; 

Fig. 5 illustrates a block diagram for a system of processing digital video with 
a subsystem that contains a computer control program according to an embodiment of the 
present invention; 

Fig. 6 illustrates a block diagram for the various functions of the computer 
control program according to an embodiment of the present invention; and 

Fig. 7 illustrates a block diagram for a system of processing digital video with 
subsystems that performs various processing steps according to an embodiment of the present 
invention. 

DESCRIPTION OF THE SPECIFIC EMBODIMENTS 

According to the present invention, techniques including methods and systems 
for digital video processing are provided. 

Fig. 1 illustrates a block diagram according to an embodiment of the present 
invention. The inventive method comprises three steps: a video providing step 10, a 
sigmaSAD calculation step 20, and a divisor calculation step 30. 

In the video providing step 10, a stream of video frames is provided. The 
video data may be derived from live video data, pre-recorded video data, or the like. Various 
video frames are assigned or labeled to various compression modes, such as I, P or B. These 
compression modes significantly affects the size of the compressed data, irrespective the 
encoding parameters that are to be selected in the next two steps. 

The SigmaSAD calculation step 20 comprises a motion estimation step 40 and 
a sum of absolute differences (SAD) averaging step 50. In the motion estimation step 40, the 
optimal motion vector for each microblock of a current frame is determined by methodically 
searching for the vector yielding the smallest difference between a current microblock and a 
microblock in a previous frame. The resulting figure of merit is called the sum of absolute 
differences (SAD). Once the SAD value for a microblock is determined, it is written to a unit 
of a reference memory. In the reference memory, each unit typically corresponds to one 
microblock, and each microblock of the current frame corresponds to one unit. For example, 



if the current frame consists of 99 quarter common interface format (QCIF) microblocks, the 
reference memory stores 99 SAD values. 

In the SAD averaging step 50, the SAD values for all of the microblocks of the 
current frame is summed up, and the sum is then divided by the total number of microblocks. 
The resulting value is called sigmaSAD. Since the sigmaSAD value of a frame must be 
calculated before the encoding of any microblocks in the frame, the first encoded microblock 
would be delayed by the duration of the motion estimation step 40. For example, if the 
encoder can process 30 frames per second, the delay would amount to only less than 20 msec. 

In the divisor calculation step 30, a divisor for each frame is determined. The 
divisor is defined as the ratio between sigmaSAD and a quantization step value (Q) or the 
ratio between sigmaSAD and the sum of Q and a constant. The constant may be set to -1 . In 
the present embodiment, SigmaSAD has been determined to be an accurate predictor for the 
frame complexity for a large number of standard QCIF sequences. Experimentally, the 
Pearson product-moment correlation between sigmaSAD and the unencoded data size of a 
frame (hereinafter "unencoded frame size") was found to be approximately 0.9. 

Even though sigmaSAD correlates with the unencoded frame size, sigmaSAD 
does not vary linearly with the unencoded frame size. Further, Q responds nonlinearly with 
the unencoded frame size. The non-linearity of sigmaSAD and that of Q closely match each 
other. More specifically, a small variation in small sigmaSAD and a small change in small Q 
can both make large changes in encoded data size, and the converse relationship holds true 
for both high sigmaSAD and high Q. Therefore sigmaSAD and Q change almost linearly 
with each other for a given encoded data size. Although the relationship between sigmaSAD 
and Q can be even more precisely described by a complicated nonlinear formula, the linearity 
or the modified linearity by replacing Q with Q plus a constant, adopted in the present 
embodiment has generated satisfactory results, where Q varies in value from 2 through 3 1 . 

The determined divisor remains constant throughout the entire frame, but it 
may change with various frames. After a divisor is selected for one frame, if a subsequent 
frame contain the same scene, i.e. the same spatial complexity, and contains changes only in 
motion complexity, the divisor usually changes very little. 

For example, the instant embodiment may employ a set of 16 potential divisor 
values, which are indexed from 0 through 1 5 and may range from a value suitable for simple 
images to one that can handle complex scenes. To encode QCIF at ten frames per second, the 
set of divisor values can be chosen from a series of quotients ranging from 64 to 264. In 
order to choose a particular divisor value, an integer is selected such that its quotient modulo 



16 is the divisor index, and this integer is called complexity integer. For every I-frame, the 
complexity integer is set to a value midrange between 0 and 255. In subsequent P-frames 
whose scenes do not change significantly, the complexity integer may be adjusted if the 
buffer fullness deviates from the targeted fullness. 

The magnitude of the adjustment increases with the difference between the 
buffer fullness and the targeted fullness. The resulting change in the complexity integer 
occurs usually gradually; hence, the divisor also changes only gradually for frames 
possessing almost the same spatial complexity. The divisor value determined through the 
complexity integer is then used to calculate the Q value for the frame. The Q value equals 
the ratio of the sigmaSAD value to the divisor value, or equals the difference between the 
ratio of the sigmaSAD value to the divisor value and a constant, where the constant may be , 
set to -1. This process would be repeated for each frame in the incoming video stream before 
each frame is encoded. 

In any event, if the buffer fullness approaches its maximum level and is 
beyond the control through the adjustment of the complexity integer, the subsequent frames 
would be "skipped," effective alleviating the fullness of the buffer. 

Fig. 2 illustrates a block diagram according to another embodiment of the 
present invention. The inventive method comprises three steps: a video providing step 1 10, a 
sigmaSAD calculation step 120, and a divisor calculation step 130. 

In the video providing step 110, the unencoded video frames are transmitted 
into an encoding system, without labeling various video frames with different compression 
modes, such as I, P or B. 

The SigmaSAD calculation step 120 comprises a motion estimation step 140 
and a SAD averaging step 150, similar in function as discussed above. During the motion 
estimation process, however, the changes of a sequence of frames in spatial complexity and 
in motion complexity can be determined with skills well known in the art. Hence the 
decision to encode a frame in the I-mode can be made before any P-frame encoding has 
actually been accomplished. This is a far more efficient processing sequence than encoding 
microblocks in a purely sequential manner. 

Fig. 3 illustrates a block diagram according to yet another embodiment of the 
present invention. The inventive method comprises three steps: a video providing step 210, a 
sigmaSAD calculation step 220, and a divisor calculation step 230. 

In the video providing step 210, a stream of video frames is provided. 
Furthermore, various video frames are assigned to various compression modes, such as I, P or 



B. These compression modes significantly affects the size of the compressed data, 
irrespective the encoding parameters that are to be selected in the next two steps. 

The SigmaSAD calculation step 220 comprises a motion estimation step 240 
and a SAD averaging step 250, substantially similar to that discussed above. 

In this embodiment, for subsequent P-frames, the complexity integer may be 
adjusted so that the size of the encoded data for various frames would have approximately the 
same size. In order to achieve approximately the same size, frames with simple spatial 
complexity and motion complexity use small Qs, which would enable the frames to possess a 
high image quality. On the contrary, complex frames use large Qs; thus the image quality is 
reduced. However, a high image quality is not important typically for complex frames 
because these complex frames typically occur after a dramatic change in scenery and when 
the human eye-brain system cannot comprehend the details of an image. 

Fig. 4 illustrates a block diagram according to yet another embodiment of the 
present invention. The inventive method comprises three steps: a video providing step 310, a 
sigmaSAD calculation step 320, and a divisor calculation step 330. 

In the video providing step 310, a stream of video frames is provided. 
Furthermore, various video frames are assigned to various compression modes, such as I, P or 
B. These compression modes significantly affects the size of the compressed data, 
irrespective the encoding parameters that are to be selected in the next two steps. 

The SigmaSAD calculation step 320 comprises a motion estimation step 340 
and a SAD averaging step 350 as discussed above. 

In this embodiment to reduce the buffer fullness, a frame immediately 
following an I-frame can be skipped, because following a dramatic scenery change in an I- 
frame, the human eye-brain system cannot comprehend images in any detail. The locations 
of the I-frames are usually determined in the video providing step 310. To compensate for 
the skipped frames, the frames immediately preceding the I-frames, instead of the I-frames 
that are coded with large Qs and possess low image quality, would be extended for one more 
frame time, where one frame time is defined as the display duration of each frame. 

In any event, if the buffer fullness approaches its maximum level and is 
beyond the control through the adjustment of the complexity integer, the subsequent frames 
would be "skipped," effective alleviating the fullness of the buffer. 

Fig. 5 illustrates a block diagram according to an embodiment of the present 
invention. The inventive encoding system 470 comprises a memory unit 480, within which a 
computer program 490 is stored. 
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The computer program 490 instructs the system 470 to perform an encoding 
process comprising three functions: a video providing function, a sigmaSAD calculation 
function, and a divisor calculation function. 

In the present embodiment, the system 470 may be a computer system, a 
custom processing chip, or other device including memory and a processor. The video 
providing function provides a stream of video frames. These frames may be from a video 
camera, from a hard disk, CD-ROM, DVD-ROM, or the like. Furthermore, the video 
providing function assigns various video frames to various compression modes, such as I, P 
orB. 

The processor may be configured to perform the sigmaSAD calculation 
function including motion estimation and SAD averaging, and the divisor calculation 
function based upon video data, as discussed above. In alternative embodiments, dedicated 
hardware may be used to provide one or more of the above calculations. 

Embodiments of the present invention include code that directs the processor 
to perform the aforementioned functions. To further illustrate the working of the computer 
program 490, an exemplary program that performs the sigmaSAD calculation function and 
the divisor calculation function is illustrated as follows: 
// Beginning of the program 
int CMp4dll::rate_control( 

int num mblk, // number of macroblocks per picture 

int BPF, // channel bit rate divided by frame rate 

int fbits, // bits used to code previous frame 

index, // buffer content, quantitized to 4 bits 

// quantization factor 
// first P-frame 

numintramblk, // number of intra MB in this picture 
sadavg) 
divindex; 
skip; 
stuff; 

// We define a bit rate threshold applied to the average bit 
// rate generated in the previous four frames. If we are below 
// the first threshold, we defer any complexity adjustment that 
// would tend to further reduce bit rate. 



int 



{ 



int 
int 
int 
int 
int 
int 
int 



*Q, 

first_p, 



// counter for scheduled MB refresh 
// index for divisor[0. .15] 
// indicator to skip a frame (buffer full) 
// indicator to stuff bits (buffer empty) 



if (first_p) // indicates I-frame has just been co 

{ 

// initialize this number at every I-frame 
if(BPF> 0x1000) 

complexity = 0x80; 

else 

complexity = 0x60; 
first_p = 0; 

else 

{ 

bits[n] = fbits; // four element array 

n++; 

n&=3; 

} 

// Conceptually we determine the Q value for the current frame 
// by the integral of the M.E, SAD (ssad) values. These values are 
// interpreted differently with different image complexities. 
// The current buffer fullness only revises our estimate of 
// the scene complexity. 

// "sad_avg" - sigma SAD, the average MB SAD value over the frame. 

// "index" = buffer content / 4096, 0 to 15, used as a selector 

// in the case statement, which modifies "complexity". 

// "complexity" - 0 to 255, a measure of the spatial data content 

// of the image being coded. With larger complexity a 

// given SAD value will generate fewer bits. Controlled 

// by buffer fullness with a long time constant. 

// Initialized after coding an I-frame. 

// "divindex" = complexity / 16, used as an index to divisor[], 
// takes values from 0 to 1 5 

// "divisor[] n - an array that defines 16 relationships between 

// image complexity and SAD integral. 64 --> 264 

// "stuff 1 - a command to insert stuffing bits 

// "skip" - a command to skip the next scheduled frame coding 
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if(index<0) 

index = 0; 
if (index > 15) 

index = 15; 
5 bitrate = 0; 

for (i=0;i<4;i++) 

bitrate += bits[i]; 

bitrate »= 2; // av. bits/frm over last 4 frames 

// We impose an immediate complexity correction when the average P-frame 
10 // bit rate over the previous four frames exceeds 10,240 bits per frame 

I // (1 .6 times the nominal rate). The limit declines linearly with respect 

W ! // to buffer fullness to 8448 bpf (1 .32 x) when the buffer is more than 

iSjj" //half full. 

(50 if(index<8) 

W limit = (96 * BPF) » 6; // 1 .5 * BPF 

Pi else // 1.5 to 1.125* BPF 

n limit = ((39 - index) * BPF * 3) » 6; 

;J if (bitrate > limit) 

111! complexity -=16; // reduces divindex by one 

|! stuff =0; 

. ,.j skip = 0; 

switch (index) 

{ 

case 0: // buffer almost empty 

25 if(bitrate<BPF) 

complexity += 64; 
stuff = 1; 
break; 

case 1: 

30 if (bitrate < BPF) 

complexity += 32; 
stuff = 1; 
break; 

case 2: 
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// Similarly, if we are above this threshold we defer any 

// complexity changes that would increase bit rate. 

int bitrate; // sum of the four values 

int limit; // sets limit on rate of filling the buffer; 

int ij, status; 

status = 0; // default - not I-frame, no stuff, no skip 

if(BPF!=oldbpf) 

{ 

j=64; 

for(i = 0; i<16;i++) 
{ 

divisorfi] = j; 
j += (BPF/divfactor); 

} 

oldbpf=BPF; 

} 

int divisortemp[16] = {64,96,112,128,136,144,152,160,168,176,184,192,200,216,232,264}; 
#ifdef DEBUG 

printf( , 3PF=%d,fbits=%d,index=%d,Q=%d,fu-st_j=%d,niun_intra_mblk=%d,sad_av 
g=%d\n", 

BPF,fbits,index,*Q,first_p,num_intra_mblk,sad_avg); 

#endif 

threshold = (num_mblk«l)/3; 
// Check if should be I picture 
if (num_intra_mblk > threshold) 
{ 

status = 1 ; 
return (status); 

} 

if(first_p) 
{ 

for (i=0;i<4;i++) 

bits[i] = BPF; 

} 
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if(bitrate<BPF) 

complexity += 24; 

break; 

case 3: 

if (bitrate < BPF) 

complexity += 16; 

break; 

case 4: 

if (bitrate < BPF) 

complexity += 8; 

break; 

case 5: 

if (bitrate < BPF) 

complexity += 4; 

if (bitrate > BPF) 

complexity -= 4; 

break; 

case 6: 

if (bitrate < BPF) 

complexity += 4; 

if (bitrate > BPF) 

complexity -= 4; 

break; 

case 7: 

if (bitrate > BPF) 

complexity -= 4; 

break; 

case 8: 

if (bitrate > BPF) 

complexity -= 8; 

break; 

case 9: 

if((*Q<31)&& (bitrate > 
complexity -= 12; 
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break; 
case 10: 

if ((*Q < 31) && (bitrate > BPF)) 
complexity — 16; 

break; 
case 1 1 : 

if ((*Q < 31) && (bitrate > BPF)) 
complexity -= 20; 

break; 
case 12: 

if ((*Q < 31) && (bitrate > BPF)) 
complexity -= 24; 

break; 
case 13: 

if ((*Q < 3 1) && (bitrate > BPF)) 
complexity -= 28; 

break; 
case 14: 

if ((*Q < 31) && (bitrate > BPF)) 

complexity -= 32; 
skip = 1 ; 
break; 

case 15: // buffer almost full 

if ((*Q < 31) && (bitrate > BPF)) 

complexity -= 64; 
skip = 1 ; 
break; 

if (complexity < 0) 

complexity = 0; 
if (complexity >= 256) 

complexity = 240; 
divindex = complexity » 4; 
*Q = sad_avg / divisortemp [divindex] + 1; 
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if (*Q<2) 

*Q = 2; 
if(*Q>31) 

*Q = 31; 

#ifdef DEBUG 

printf(" Chose Q=%d, complexity=%d,divindex=%d,recent avg bits/frame=%d\n", 
*Q,complexity,divindex,bitrate); 

#endif 

if (stuff) 

status = 2; 
if (skip) 

status = 3; 

// printf ("index %2d cmplx %3d sad %5d Q %2d stat 
%d\n" ,index,complexity ,sad_avg, * Q.status) ; 
return (status); 

} 

// End of the program 

It is understood that the examples and embodiments described herein are for 
illustrative purposes only and that various modifications or changes in light thereof will be 
suggested to persons skilled in the art and are to be included within the spirit and purview of 
this application and scope of the appended claims. All publications, patents, and patent 
applications cited herein are hereby incorporated by reference for all purposes in their 
entirety. 
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